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Green's functions are a central element in the attempt to understand non-perturbative phenomena 
in Yang-Mills theory. Besides the propagators, 3-point Green's functions play a significant role, 
since they permit access to the running coupling constant and are an important input in functional 
methods. Here we present numerical results for the two non-vanishing 3-point Green's functions in 
3d pure SU{2) Yang-Mills theory in (minimal) Landau gauge, i.e. the three-gluon vertex and the 
ghost-gluon vertex, considering various kinematical regimes. In this exploratory investigation the 
lattice volumes are limited to 20 3 and 30 3 at j3 = 4.2 and j3 = 6.0. We also present results for 
the gluon and the ghost propagators, as well as for the eigenvalue spectrum of the Faddeev-Popov 
operator. Finally, we compare two different numerical methods for the evaluation of the inverse of 
the Faddeev-Popov matrix, the point-source and the plane-wave-source methods. 
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I. INTRODUCTION 

The non-perturbative properties of Yang-Mills theory 
are still an open and challenging problem, especially the 
issue of confinement. Nonetheless, much progress has 
been made in their understanding over the decades. One 
central element in these investigations are Green's func- 
tions, which can describe a quantum field theory com- 
pletely. In particular, their infrared behavior has been 
related to confinement: in two of the most popular sce- 
narios of confinement, the Gribov-Zwanziger scenario 
tl and the Kugo-Ojima scenario a specific 
behavior is predicted for the 2-point Green's functions 
— the propagators — in Landau gauge. 

According to these scenarios, the Faddeev-Popov ghost 
propagator should be enhanced in the infrared limit when 
compared to a massless-particle pole, while the gluon 
propagator should vanish. Calculations using functional 
methods agree with these predictions and have found 
that the propagators show a power-like behavior in the 
far infrared region, with characteristic exponents. These 
results have been obtained usingDyson-Schwinger equa- 
tions in various dimensions 0,0, 111 and have been 
confirmed by renormalization-group methods . On the 
other hand, these methods rely on approximations con- 
cerning the higher n-point Green's functions, especially 
the three-point vertices. Thus, it is important to check 
explicitly whether the assumptions made for the vertices 
in these methods are justified. This has been done, again 
using functional methods, but as of yet only for specific 
momentum configurations of the n-point Green's func- 
tions and in the far infrared limit |lCt 1 1 1 1] , or using fur- 
ther approximations [Til ll2T| . In both cases the results 
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are consistent with the assumptions made so far. 

Numerical studies of lattice gauge theories also support 
the two confinement scenarios described above, confirm- 
ing the infrared enhancement of the ghost propagator 

[Hill EE E3 

and the suppression of the gluon prop- 
agator (in the 3d case) at low momenta |lj], Ujj- Let 
us recall that, in order to probe the infrared limit, one 
needs to use very large lattice volumes. Recently, in four 
dimensions, investigations have been performed also con- 
sidering (strongly) asymmetric lattices [Tflj . in order to 
have access to small momenta while keeping the lattice 
volume relatively small. The results obtained show an 
infrared suppression of the gluon propagator. However, 
these simulations are affected by systematic effects |20| . 
making it difficult to extract quantitative information 
from this type of lattices. 

Let us stress that the agreement between the functional 
methods and the lattice data is still at the qualitative 
level. In particular, it is still not clear if the numer- 
ical data su ppo rt a gluon propagator vanishing in the 
infrared [TsL l2lj|. Recent studies by functional methods 
[22| suggest nontrivial effects related to the use of dis- 
cretized space-time on compact manifolds, which could 
be responsible for this discrepancy. Let us also recall 
that, in the continuum, a finite (non-zero) gluon propa- 
gator at zero momentum seems to be compatible with the 
Gribov-Zwanziger scenario only in the case of an infrared 
enhanced ghost-gluon vertex |23i l24| . On the other hand, 
such an enhancement would be at variance with results 
from functional methods |ToL ITTl IT2] ] and from recent lat- 
tice studies [2fJ, |26j . 

In this work we present an exploratory study of the 
Landau gauge three-gluon vertex and of the ghost- 
gluon vertex for various momentum configurations, in 
the three-dimensional pure SU(2) case. After this ex- 
ploratory study we will be able to consider (for the in- 
teresting cases) very large lattice volumes and to probe 
the limit of small momenta. This is the first numeri- 
cal study of these vertices in the three-dimensional case. 
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Clearly, the 3d simulations are computationally less de- 
manding than those in 4d and they may help to get a 
better understanding of the four-dimensional case. Fur- 
thermore, results in three dimensions are of interest in 
their own right, as three-dimensional Yang-Mills theory 
is the (most relevant part of the) infinite-temperature 
limit of its four-dimensional counterpart |27|. In addi- 
tion, a connection exists between 4d Yang-Mills theory 
in Coulo mb g auge and 3d Yang-Mills theory in Landau 
gauge [UHg. 

Let us note that previous lattice studies of these ver- 
tices in the 4d case j2^, HjJ, 1^3, usually focused on 
specific momentum configurations, mostly with the aim 
of extracting the running coupling constant and for com- 
parison to perturbative studies. In this preliminary work 
we use several kinematical configurations in order to test 
the numerical methods employed and to study the in- 
fluence due to discretization and to finite- volume effects. 
In addition, we present results for the propagators and 
for the spectrum of the Faddeev-Popov operator, which 
plays an important role in the Gribov-Zwanziger scenario 

Eillllii. 

In Section [H] the technical details for the generation 
of the configurations, the gauge-fixing procedure and the 
error analysis are given. We report definitions and results 
for the propagators in Section IlIII and for the vertices in 
Section llVI A summary and outlook conclude this work 
m Section |V| 



II. GENERATION OF CONFIGURATIONS 

The action considered is the usual, unimproved Wilson 
action [3i| for the SU(2) gauge group. Configurations are 
generated using a hybrid-over-relaxation (HOR) update, 
consisting of five over-relaxation 34] and one heat-bath 
sweeps. For the heat-bath update, a mixed Creutz [3{| 
and Kennedy-Pendleton j3(| algorithm is employed. The 
lattices have volumes V — N 3 = 20 3 and 30 3 , and cal- 
culations have been performed at j3 = 4.2 and = 6.0. 
We use 200 HOR updates for thermalization and 40 or 
45 HOR updates between evaluation of the Green's func- 
tions. The results have been obtained in multiple inde- 
pendent runs using hot initial configurations. 

For the extraction of the vertices, large statistics have 
been necessary. Table[|]lists the precise values. Also, the 
expectation value of the plaquette is given for each lattice 
volume and (3 value. The results are in agreement with 
[l8j and [37j . The error on the plaquette is the statistical 
error including a correlation-time analysis. We always 
find that the integrated auto-correlation time |3^| is less 
than 1 HOR sweep. Following Ref. we a l so evaluate 
the inverse lattice spacing (in GeV). 

The gauge fixing to Landau gauge has been performed 
using a stochastic-over-relaxation method 39] with a self- 
adapting acceptance probability. The condition for gauge 



TABLE I: For each lattice volume V and coupling j3 we report 
here the average value of the plaquette < P >, the inverse 
lattice spacing 1/a (in GeV), the number of independent runs 
and the number of configurations considered for the evalua- 
tion of propagators and vertices. The values of the plaque- 
tte, the data for the propagators and the eigenvalue spectrum 
(reported in Section IllH have been obtained using the data 
set labeled Propagators. On the other hand, the propagators 
involved in the determination of the vertices have been de- 
termined considering the same data set as the corresponding 
vertex. Finally, note that the set labeled Three-gluon ver- 
tex includes all the configurations considered in the other two 
data sets. 
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fixing has been a test on the quantity 

-y-y 

A AT 



[tr(Q^ c )] 2 
£(tr{fe.(*/i) ~ Q>J) 2 , 



(1) 



which was required to be less than 10 12 . Here we use 
the definitions 



W = Yi y [9(x)Up(x)g(x + erf 



-g(x + e^)U u (x^g(x^] (2) 



(3) 



Also, {U^(x)} is a thermalized lattice configuration, 
{g (x)} represents the gauge transformation applied on 
the link variables U^x), the symbol t indicates Hermi- 
tian conjugate, is the lattice side in the \i direction, d 
is the space-time dimension (3 in our case), e M is a pos- 
itive unit vector in the fi direction and a c are the three 
Pauli matrices, normalized as c 2 = 1. Both the link 
variables U^x) and the gauge-transformation matrices 
g(x) are elements of the SU(N C ) group (in the funda- 
mental iV c x N c representation). In order to speed up 
the calculation, the quantity ee was not evaluated for 
each stochastic-over-relaxation sweep, but an adaptive- 
predictor method has been used for each configuration. 
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Thus, in some cases we obtain for eg a value consider- 
ably smaller (by some orders of magnitude) than 10 -12 . 
Nevertheless, we usually found for the quantity e§ a final 
value between 10~ 12 and 1CP 13 . No results on Gribov- 
copy effects will be given in this exploratory study, but 
tests using relatively small statistics, i.e. generating 20 
Gribov copies for each configuration and considering only 
a small subset of thermalized configurations, suggest that 
the effects are at the quantitative rather than at the qual- 
itative level. 

All errors given for the propagators and the vertices 
have been calculated using a standard bootstrap method 
with 1000 bootstrap samples. The quoted errors repre- 
sent a 67% confidence interval. 

Most of the results shown here have been obtained us- 
ing a code developed independently from the one used in 
Ref. (see also Nonetheless, when possible, we 

checked that results obtained with the two codes agree 
for propagators and vertices. 



III. PROPAGATORS 

A. Gluon propagator 

The gluon propagator is given by the correlation func- 
tion (see for example [lfij ) 

d;1(p) = ^<ai(p)aK-p)>, (4) 

with the momentum-space lattice gluon field defined as 

m 

X 

(5) 

Here the components p^ of p get the integer values 
—N/2 + 1 , . . . , N /2 (for even lattice sides and for sym- 
metric lattices with = N) . Note that we do not divide 
the Fourier transform of the gluon field by the lattice vol- 
ume V = N 3 . 

As for the exponential pre-factor exp^iirp^/N) = 
exp(i7rppa/L), it allows one to obtain an improved lattice 
Landau gauge condition 30] , i.e. the continuum condition 
d ■ A = is recovered with corrections 0(a 2 ), instead of 
the corrections O(a) obtained when this pre-factor is ne- 
glected. This exponent also appears naturally in a weak- 
coupling expansion of lattice gauge theory |4j . Note that 
this factor is a discretization correction. Indeed, it goes 
to 1 in the limit a — > while keeping the physical lattice 
size L = aN fixed. Furthermore, it is equal to 1 in the 
infrared limit p — ► 0. This factor cancels when evaluat- 
ing the scalar part of the gluon propagator [see Eq. © 
below], but in general not for vertices. An exception is 
e.g. the orthogonal configuration of the ghost-gluon ver- 
tex considered below (see Section llVll . Numerical studies 
[3fl E2I] have verified that this factor is also necessary in 



order to obtain the correct tensor structure of Green's 
functions involving gluon fields. 

After contracting Eq. with a transverse projector 
and a unit matrix in color space, the scalar part of the 
gluon propagator is given by 

D{ p) = vjfT,< [^»] 2 + [* A »] 2 >, ( 6 ) 

where 5L4° (p) and SL4^(p) indicate, respectively, the real 
and the imaginary part of A®(p) and the normalization 
Af is given by dN c for p > and by (d — 1)N C for p = 0. 
The propagator is by definition inherently positive semi- 
definite. Let us recall that in minimal Landau gauge one 
has 

EW)=0, (?) 

where the components of the physical momenta (denoted 
by capital letters when they are in lattice units) are given 
by 

P„ = 2 S in^. (8) 

Results will be presented as a function of the magni- 
tude of the physical momentum p = |P|/a (in GeV). 
Note that, with our notation, the continuum gluon prop- 
agator is obtained considering the product (3a 2 D(k). 
Indeed in d dimensions one has (3 = 2N c /(g 2 a 4: ~ d ). 
Also, with N c = 2, the lattice quantity 2A®(x)/(ga) = 

y/ (3 /a d ~ 2 A'f l (x) goes to the continuum quantity A^(x) 
in the formal continuum limit a — > 0. In the same limit, 
\J (3 a d+2 j4°(p) converges to the continuum momentum- 
space gluon field A"(p). Thus, for any dimension d, 
the lattice quantity f3a 2 D{k) converges to the continuum 
gluon propagator in momentum space. 

It is necessary to evaluate the gluon propagator for 
various momentum configurations in order to determine 
the vertices. In addition, some extra momentum con- 
figurations have been considered to check effects related 
to the breaking of rotational invariance. The various 
kinematical configurations are given in Table [H] Let 
us note that, in general, for a given momentum p one 
has A^ljt) 7^ A"(p), when \i ^ v. For example, if the 
momentum is aligned along the /i-direction, then A"(p) 
vanishes for each color component due to the relation 
(0) or, equivalently, due to the constancy of the quantity 
q^Xfj) defined in Eq. 

The results for D{p)p 2 and D(p) are shown in Fig. 
^ for the 30 3 lattice volume and the two (3 values con- 
sidered. The statistical errors are below 1%, thus sys- 
tematic errors are clearly visible. In particular, violation 
of rotational symmetry is the dominant effect and is of 
the order of several percent. In addition, finite-physical- 
volume effects are visible, especially in the far infrared, 
when going from (3 = 4.2 to (3 = 6.0. The latter effect is 
much more pronounced for the 20 3 lattices. Altogether, 
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TABLE II: Kinematical momentum configurations consid- 
ered for the propagators. The columns Gluon and Ghost iden- 
tify whether these configurations have been used for the eval- 
uation of the propagators. The columns n x , n y and n z give 
the components of the momentum. The Type is used for later 
identification of classes of momenta. The quantities p m i n and 
Pmax (in MeV) are the smallest and the largest momenta. The 
smallest momentum is evaluated at f3 = 4.2, while the largest 
is evaluated at (3 — 6.0 (in both cases we used the 30 3 lattice 
volume). The variables n and m run independently over all 
possible positive integer values 0, . . . , N/2 (but for the ghost 
propagator one cannot consider the zero momentum case). 
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the systematic errors are the dominating effect at this 
level of statistics. 



At the largest momentum (for (3 — 6.0) the propaga- 
tor is about 6% or 5.5% larger than its re-summed one- 
loop perturbative value [8j (which itself deviates from 
the tree-level value by about 3.5%) on the 20 3 and the 
30 3 lattices, respectively. This discrepancy cannot be ex- 
plained by removing tadpole contributions from the def- 
inition of the gluon field (see for example ) . Indeed, 
if we consider ^3] the gauge invariant definition of the 
tadpole factor given by u =< P > 1 / 4 , where < P > is 
the average value for the plaquette, we find (see Table UJ) 
that the gluon propagator is multiplied by ~ 1-1610 
at j3 = 4.2 and by w 1.1011 at (3 = 6.0. This 

clearly enhances the discrepancy between the data and 
the re-summed leading-order perturbation theory. On 
the other hand, one should recall that, by changing the 
lattice discretization for the gluon field (see for example 
|44j and references therein), one finds a gluon propaga- 
tor that differs by a global multiplicative constant. In 
particular, one can easily find definitions of the gluon 
field for which the gluon propagator is sensibly smaller 
than that obtained with the standard definition (see for 
example Table 1 in Ref. 0]). Of course, different lat- 
tice discretizations converge to a common result as the 
continuum limit is approached |45| . 

Finally, we note that the combinations of lattice size 
and (3 values considered here are not sufficient to reach 
the infrared regime, where the bending over of the prop- 
agator has been observed Nevertheless, one clearly 
sees from Fig.^that the propagator is less singular than 
1/p 2 for momenta smaller than about 1 GeV. 
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FIG. 1: The gluon dressing function D{p)p 2 (top) and the 
gluon propagator D(p) (in GeV -2 , bottom) as a function of 
p (in GeV). Open and full symbols correspond respectively to 
f3 = 4.2 and to (3 = 6.0. In both cases we consider the lattice 
volume V — 30 3 . Squares and triangles correspond to plane 
and diagonal momentum configurations, respectively. 



B. Ghost propagator 

The numerical evaluation of the ghost propagator is 
considerably more complicated than that of the gluon 
propagator. Indeed, one has to evaluate 



Dg(p) = ~ < (M~ 1 ) ab (p) >, 



(9) 



where M ab (x,y) is the Faddeev-Popov operator, defined 
in the continuum as 



-d^Df = S(x - y)(-d 2 6 ab + gf^d^j 



(10) 



Here f abc are the structure constants of the SU(N C ) 
gauge group. In Landau gauge, on the lattice, this op- 
erator is a matrix (with color and space-time indices), 
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defined by its action on a scalar function uj h (x) (with 
color index b) as 2] 

M ah ( x , y)u>\y) = s xy Y^{Gt(y)l" h (y) - ^ b (y + ^)] 

c 

-^(y-^K^-e^)]}. (11) 

Here, the sum over repeated indices (y and b) is under- 
stood and G ab (x) is given by 

Gf(x) = \tr{{a a ,a b ) [U^x) + U^{x) + ]) , (12) 

i.e. it is proportional to S ab . Let us recall that in minimal 
Landau gauge 2\ the gauge-fixed configurations are in- 
side the first Gribov horizon, implying that the Faddeev- 
Popov operator M ab (x,y) is symmetric and positive (in 
the subspace orthogonal to the trivial and constant zero 
modes). 

In order to evaluate the Fourier transform of the in- 
verse operator 1 

(M~ 1 ) ab (p, q) = J2 e2m{pX+9V)/N ( A ' rl ) ab ( x >y) ( 13 ) 
x,y 

the matrix inversion has been performed using the point 
source S ac (5 X Q — (see Ref. |46|). Compared to the 
inversion using a plane- wave source |47j , this method has 
the advantage of only N% — 1 inversions per configura- 
tion, independently of the number of momenta, instead of 
N% — 1 inversions per configuration and for each momen- 
tum considered. However, as we will see below, statistical 
fluctuations are significantly enhanced at large momenta 
|46j. Also, note that this procedure is ambiguous with 
respect to the sign of the resulting propagator 2 (or of 
the eigenvalues). Thus, the sign has to be assigned by 
hand. 

Since the Faddeev-Popov operator M ab (x,y) is sym- 
metric and positive, the matrix inversion can be per- 
formed u sing a conjugate-gradient (CG) method (see for 
example |48j|). From the numerical point of view, there 



1 In the general case one considers p ^ q. However, in the case 
of the ghost propagator, one has p = —q due to momentum 
conservation. Note that, again, we do not divide by the lattice 
volume V when considering the Fourier transform. 

2 When considering the plane-wave source the ghost propagator 
is the expectation value of a positive operator on the plane- 
wave state, i.e. the result is always positive-definite. This is 
not true with the point-source method. Indeed, one can obtain 
(on a given lattice configuration) a negative value for the quan- 
tity (M~ 1 ) ab (p). Nevertheless, on average, the ghost propagator 
Dq(p) is positive for all values of the momentum p also when us- 
ing the point-source method. 



are two issues to be careful about. First of all, as finite- 
precision arithmetic is used, it is possible for the solution 
to develop a component along the (constant) zero modes 
of the Faddeev-Popov operator. Thus, it is necessary to 
reorthogonalize the solution at each CG iteration with 
respect to this subspace. 

The second aspect is related to the convergence of the 
CG method. Indeed, in finite precision arithmetic, the 
CG method normally loses orthogonality to the Krylov 
subspace already spanned. Thus, the magnitude of the 
residual can be incorrect. To compensate for this, an ad- 
ditional, albeit technical, quality criterion can be consid- 
ered (see ji^l section 5.3 for details). As a convergence 
test, we check if the average of the components of the 
residual is less than 10~ 12 and if the additional quality 
parameter 3 is less than 1CP 13 between two CG itera- 
tions. This is usually achieved after several tens to a few 
hundreds of iterations. We have also checked that much 
weaker convergence tests do not change the propagators 
by more than a few percent. Let us note that one CG in- 
version of the Faddeev-Popov matrix is much faster than 
the gauge-fixing procedure when considering large lattice 
volumes. 

The results for the color-averaged diagonal ghost prop- 
agator are shown in Fig.[2]for the lattice volume V = 30 3 . 
Again, the statistical error is essentially negligible com- 
pared to the systematic effects. Note that the conse- 
quences of violation of rotational symmetry are much less 
severe than in the gluon case, although they are still the 
dominant source of systematic errors. As the ghost is 
a scalar, this was somehow expected. Also, the finite- 
physical- volume errors are somewhat smaller than in the 
gluon case. Comparing the lattice data to the perturba- 
tive predictions y| at the largest momentum (for (3 = 6.0) 
we find that the former are larger by about 6% and 5.5% 
for the lattice volumes V = 20 3 and V — 30 3 , respec- 
tively. In this case, perturbation theory only deviates 
about 1% from the asymptotic tree- level value at this 
momentum. If one considers tadpole improvement |15| . 
then the ghost propagator gets multiplied by uq, i.e. by 
0.92807 at j3 = 4.2 and by 0.95298 at /3 = 6.0, leading to 
a better agreement with perturbation theory. 

Here we didn't try to fit the data for the ghost prop- 
agator and obtain an estimate for its infrared exponent. 
Nevertheless, from Fig. [21 one clearly sees that Dc{p) is 
more singular than 1/p 2 . 

We also compared the results obtained using the two 
different inversion methods mentioned above (i.e. point- 
source and plane- wave-source methods). In the second 
case the inversion has been done using a CG method with 
even-odd preconditioning and a convergence test given by 
l r | 2 /l r o| 2 < 10 -8 , where r and ro are the final and the ini- 
tial residuals, respectively. The use of the even-odd pre- 



In the present case, due to memory restrictions, we only consider 
subsequences of the CG iteration of length 1. 
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Ratio of the ghost propagators from point and plane-wave source 
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FIG. 2: The color- averaged diagonal part of the ghost dress- 
ing function Dc(p)p 2 (top) and of the ghost propagator Dg(p) 
(in GeV -2 , bottom) as a function of the momentum p (in 
GeV). Open and full symbols are used for /3 = 4.2 and /3 = 6.0, 
respectively. In both cases we consider the lattice volume 
V = 30 3 . Squares and triangles correspond to plane and di- 
agonal momentum configurations, respectively. 



FIG. 3: In the top panel we show the ghost dressing func- 
tion Dg{p)p 2 obtained using the plane- wave source (full sym- 
bols) and the point source (empty symbols) for V = 30 3 and 
j3 = 4.2. Note the logarithmic scale on the y axis. In the bot- 
tom panel we report the ratio of the point-source data (PS) 
with the the plane-wave-source results (PW) for the ghost 
propagator. In both panels, circles denote momenta of type 
(0, 0,p), triangles indicate diagonal momenta and the quanti- 
ties are considered as a function of the momentum p (in GeV) . 
Here, for both methods, we considered 380 configurations. 



conditioning usually reduces the number of CG iterations 
by about a factor 2. We find that the results from these 
two methods agree on the average for the ghost propaga- 
tor as well as for the ghost-gluon vertex, discussed below. 
As can be seen in Fig. |3 the ghost propagator obtained 
using the point source oscillates around the (smoother) 
result obtained using the plane- wave source. This oscilla- 
tory behavior is stronger for the ghost-gluon vertex, but 
vanishes in both cases when the statistics is increased, 
albeit much more slowly for the latter case. 

An important assumption in functional calculations is 



that the ghost propagator is color diagonal 4 . In Fig. 
0] we show results for the real part of the off-diagonal 
components of the ghost propagator 5 . The majority 



4 We also checked that the real part of the off-diagonal components 
of the gluon propagator is essentially zero within statistical errors 
and that its central value decreases with increasing statistics. Of 
course, the imaginary part vanishes identically. 

5 Let us recall that the fluctuations of the imaginary part of the 
off-diagonal elements of the ghost propagator are connected to 
the possible existence of a ghost condensate |4!.I . In particu- 
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FIG. 4: The real part of the color upper-triangular matrix 
elements of the ghost dressing function Dq" (p)p 2 as a function 
of the momentum p (in GeV). Open and full symbols indicate 
j3 = 4.2 and /3 = 6.0, respectively. In both cases we consider 
the lattice volume V = 30 3 . Squares and triangles correspond 
to plane and diagonal momentum configurations, respectively. 



FIG. 5: A histogram of the real part of the off-diagonal 
elements of the ghost propagator (in lattice units) evaluated 
at the smallest momentum for the lattice volume V = 30 3 at 
/3 = 4.2. We also plot a Gaussian with the same mean and 
standard deviation as the histogram of the data. 



C. Eigenvalue spectrum of the Faddeev- Popov 
operator 



of the points are compatible with a null result, within 
the 67%-confidence interval, and the mean value of the 
(real part of the) off-diagonal propagator decreases for 
all momenta with increasing statistics. On the other 
hand, these fluctuations exhibit significantly enhanced 
tails (see Fig. EJ. We checked that these large fluctu- 
ations are in most cases related to large values of the 
color-diagonal part of the ghost propagator and to con- 
figurations for which gauge fixing to Landau gauge re- 
quired many more iterations than in the average case. 
Thus, these tails could be related to the exceptional con- 
figurations observed in Il6[. Let us note that, by using 
the point-source method 46] for the inversion, one needs 
very large statistics in order to see a null average for 
the off-diagonal points in the infrared. Also, using this 
method, the real (respectively, imaginary) part of the in- 
verse Faddeev-Popov matrix (M~ 1 ) ab (p, —p) is symmet- 
ric (respectively anti-symmetric) only on average and not 
for each lattice configuration, as is the case when using 
the plane- wave source |49( . 



The CG method has a close relationship with the so- 
called Lanczos algorithm, which can be used to extract 
the eigenvalues of a matrix 48]. As a consequence, it is 
possible to determine the ghost propagator and simulta- 
neously obtain information on the eigenvalue spectrum. 
Note that only the exact CG algorithm is guaranteed to 
obtain the correct spectrum (up to degeneracy) of the 
Faddeev-Popov matrix after at most (JV| — 1)V inver- 
sion steps. On the other hand, working in finite-precision 
arithmetic (see Section 4 in 0|)> the algorithm typi- 
cally does not deliver a sufficiently accurate numerical 
approximation to all eigenvalues of this matrix. Indeed, 
as the number of iterations increases, the method has 
the property that the extremal eigenvalues evaluated are 
progressively improved approximations of the extremal 
eigenvalues of the matrix considered |5l). At the same 
time, since the procedure is stopped before all eigenval- 
ues are found, the middle of the spectrum is usually not 
reliable and under-populated. Moreover, each eigenvalue 
will be found only once, independently of its degeneracy. 
Still, when averaging over many configurations — i.e. 
when considering a histogram of the eigenvalues found 
for all the configurations, normalized by the total num- 
ber of eigenvalues 6 — one should be able to obtain an 



lar, non-Gaussian fluctuations could indicate the existence of a 
spontaneous-symmetry breaking and a non-vanishing value for 
this condensate. (For a different interpretation, see Ref. HJ.) In 
this work we do not present data for the ghost condensate. 



Doing this normalization for each inversion and each configura- 
tion separately does not yield a different result within the statis- 
tical errors. 
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idea of the density of the eigenvalues at both ends of 
the spectrum. In particular, one can verify whether the 
shape of the eigenvalue density (at small eigenvalues) is 
flat or increases with the magnitude of the eigenvalues u>, 
thus gaining information on the validity of the Gribov- 
Zwanziger scenario. On the other hand, the upper end 
of the spectrum, which is also obtained, is essentially de- 
termined by perturbative contributions and thus not as 
interesting. 



Full spectrum of the Faddeev-Popov operator 
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x10' 3 




to 



FIG. 6: The density p{u) of the eigenvalues of the Faddeev- 
Popov operator (using the Lanczos method) as a function of 
their magnitude u) (in lattice units), normalized by the total 
number of eigenvalues. In the top panel we show the complete 
results while in the bottom panel only the density for small 
eigenvalues is presented. In both cases we considered 100 bins 
(this explains the different scale on the y axis). The solid 
line (in both panels) refers to the lattice volume V = 30 3 at 
P — 4.2. In the bottom panel, the dashed line correspond to 
V = 30 3 at (3 = 6.0, the dotted line to V = 20 3 at p = 4.2 
and the dashed-dotted line to V = 20 3 at f3 = 6.0. The 
total number of eigenvalues determined are 4374822 (V = 20 3 , 
P = 4.2), 11312874 (V = 30 3 , P = 4.2), 3549202 (V = 20 3 , 
P = 6.0) and 9862721 (V = 30 3 , P = 6.0). 



TABLE III: Largest and smallest eigenvalues of the Faddeev- 
Popov matrix (respectively, u>i and lu s , both in lattice units) 
for each lattice volume V and coupling p. We also report 
the corresponding eigenvalues of the lattice Laplacian (for the 
same lattice side N). For the lattice volume V = 20 3 we 
consider 200 configurations at P — 4.2 and 400 at /3 = 6.0; for 
V — 30 3 we have 380 and 500 configurations at /3 = 4.2 and 
6.0, respectively. 



Volume 


P 








Us.Lapl. 


20 3 


4.2 


11.1218(9) 


0.0072(2) 


11.9261 


0.0246 


30 3 


4.2 


11.1337(5) 


0.00275(6) 


11.9671 


0.0110 


20 3 


6.0 


11.4050(4) 


0.0102(2) 


11.9261 


0.0246 


30 3 


6.0 


11.4149(2) 


0.00353(7) 


11.9671 


0.0110 



| Smallest eigenvalue of the Faddeev-Popov operator | 




1/L 



FIG. 7: Plot of the smallest eigenvalue u s (in GeV 2 ) of the 
Faddeev-Popov matrix, as a function of the inverse lattice side 
1/L (in GeV), and of the fitting function a/L c with a — 18.127 
and c = 2.6152. 

The results are shown in Fig. (top panel). The full 
spectrum is shown for completeness. As said above, the 
small density of the spectrum for intermediate magnitude 
u> of the eigenvalues is very likely an artifact of the algo- 
rithm. For very small magnitude u> (see Fig. bottom 
panel) the spectrum shows a linear increase with lo. The 
steepness is much larger on the larger lattices and it also 
seems to increase with the physical volume, i.e. as (3 de- 
creases. These results are in qualitative agreement with 
those reported in Ref. [3i| for the four-dimensional case. 
Note, however, that a linear increase at small momenta 
is different from the result obtained in Coulomb gauge 
in four dimensions |52| , where a power- like behavior has 
been observed. 

Using a preconditioned CG method we have also eval- 
uated the smallest and the largest eigenvalues of the 
Faddeev-Popov matrix M . Results are reported in Table 
IIIII For comparison we also give the smallest and the 
largest eigenvalues of the lattice Laplacian —A for the 
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same lattice size. Let us recall that in d dimensions, the 
smallest eigenvalue of —A is given by 4sin 2 (7r/2iV), while 
the largest is equal to 4<isin 2 [7r(iV — 1)/2N}. We see that 
the largest eigenvalue of M depends weakly on the lattice 
volume V and on the the coupling (3. On the other hand, 
for a given (3, the smallest eigenvalue decreases faster 
than the corresponding eigenvalue of —A, as the lattice 
side increases. This is in agreement with the results re- 
ported in ■ Also, if we consider the smallest eigenvalue 
of M in physical units, i.e. if we multiply it by 1/a 2 , we 
find the results reported in Fig.[7| If one tries a fit to data 
using the Ansatz a/L c we obtain that the data are well 
fitted for c = 2.62(8) (with x/d.o.f. = 1.4). This seems to 
suggest that for these lattice volumes and (3 values we are 
in the scaling region for to s and that this eigenvalue goes 
to zero when the infinite- volume limit is approached. As 
a consequence, in the continuum limit, the average lattice 
Landau configuration should belong to the first Gribov 
horizon, supporting the Gribov-Zwanziger mechanism of 
confinement. 



IV. VERTICES 

In Yang- Mills theory (without valence and sea quarks) 
there are two non-vanishing 3-point Green's functions: 
the three-gluon vertex and the ghost-gluon vertex. Com- 
pared to the propagators, these are quite complicated 
since there are now three kinematical variables. In addi- 
tion, for the three-gluon vertex there exist a large num- 
ber of tensor structures [53j. Therefore, it is impractical 
to investigate all tensor structures for all possible kine- 
matical configurations. Furthermore, due to kinematics, 
the interesting signal — i.e. an appropriately defined di- 
mensionless function — is usually suppressed at large 
momenta as 1/p 6 instead of as 1/p 2 (except for some 
particular momentum configurations for which it is only 
suppressed by 1/p 4 ). Thus, it is in general quite difficult 
to obtain a good signal/noise ratio at large p. 

Let us recall that, due to momentum conservation, the 
three momenta of the vertices always lie in one plane and 
are completely characterized by the size of two of the mo- 
menta and by the angle </> between them. Two particular 
momentum configurations are of special importance. In 
the first case two of the momenta are orthogonal with 
respect to each other (this configuration will be called 
here the orthogonal one). This configuration enters into 
loop integrals with maximum angular weight 7 and, as 
we will see below, has the advantage of relatively small 

7 Let us recall that this angular weight is proportional to sin(</<) in 
3d. 



statistical fluctuations. In this case we take the two or- 
thogonal momenta respectively along the x and y axis, 
considering all possible magnitudes of both momenta in- 
dependently. If one of these two momenta vanishes, then 
one recovers the kinematical configuration used in several 
previous lattice studies of these vertices [2j| |2j| |3(], . 

In the second configuration (called here the equal one) 
the three momenta have the same magnitude and the 
angle <fi is equal to 7r/3. The vertex is then a function 
of only one variable. For this configuration the infrared 
behavior of the n-point Green's functions can be taken 
to be the limit where all n momenta vanish in the same 
way, i.e. only one scale has to be considered. In this 
case there are predictions for the behavior of the n-point 
Green's functions in the infrared limit, from studies using 
Dyson-Schwinger equations (DSEs) [l(|, [TlJ . Of course, 
on the lattice, it is not possible (in general) to select three 
equal momenta, for example, in the x—y plane. However, 
one can consider the momenta (p, — p, 0), (— p, 0,p) and 
(0,p, —p), which all have the (lattice) size V2P with P = 
2sin(7rp/iV). 

Let us recall that, due to the invariance of the lat- 
tice theory under hyper-cubic transformations, the re- 
sults should be invariant by reflection of the momenta. 
This implies that — with the exception of p = and 
of the mid-point on even lattices — we can average over 
two different kinematical configurations for each physical 
momentum. For the vertices, due to momentum conser- 
vation, we can apply the reflection transformation only 
to the independent momenta. In particular, in the or- 
thogonal case, we can apply the reflection independently 
on the two orthogonal momenta, i.e. — for each physi- 
cal momentum configuration — we can average over four 
different kinematical configurations. On the contrary, in 
the equal-momenta case, the reflection has to be applied 
at the same time on the three momenta, i.e. each time 
we can average over two different kinematical configura- 
tions. Let us note that this averaging allows us to cancel 
exactly, on each lattice configuration, contributions that 
would otherwise be zero only on average, yielding purely 
imaginary vertices. 

Finally, one has to consider the tensor structure of the 
vertices, which is not straightforward for the three-gluon 
vertex. In functional-method studies it is particular in- 
teresting to know how different a vertex is from its tree- 
level structure. Thus, we consider here the projection of 
the vertices on their tree-level values, obtaining a scalar 
function normalized to 1 when the tree-level vertex and 
the full vertex are equal. 
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A. Three-gluon vertex 

The three-gluon vertex r* 1 ^ ' abc (p, q, k) (with k = —p — q) at tree-level in the continuum has the form 

r^f' afcc (P, q, k) = ~igf abc [(q - k)n6 vp + (A - p) u S w + ( P - q) p 5^} , (14) 
while on the lattice one finds 
r%' A3 > abc ( P , q, k) = -igf abc e **d*+«»+*p)/" [(pfc)^ CO s(p,) + {k^p) v b pp cos(q p ) + (f^q) P 8^ cos(fc M )] . (15) 

Note that the cosine factors in the above equation go to 1 in the formal continuum limit a — > 0, for a fixed physical 
lattice side L — Na. The same applies to the exponential pre-factor. Here we used the notation 

p p = 2sin(p jU ) = P p (16) 

= if (17) 

with p p taking values , . . . , N/2 . Clearly, the three-gluon vertex is totally symmetric under the simultaneous 
exchange of color index, Lorcntz index and momentum, and thus completely Bose-symmetric. 

The full three-gluon vertex ^ A l/p abc is not directly available on the lattice, but one can evaluate the corresponding 
full Green's function |29| 

iA 3 ,abcf m „ j„\ _ _j_ ^ \a;^ \ b i 

V 

As said above, momentum conservation requires k = —p — q. Due to the vanishing of any vector condensates in 
Yang-Mills theory, this Green's function equals the connected Green's function, but it is still necessary to amputate 
it. The relation of the full vertex r^.^ 6 ^ (p, q, k) with the Green's function is then given by 

G^f c {p, Q, k) = Dfi(p)Dl%(q)D%(k)rfj ef (p, q, k). (19) 

In Landau gauge, one can extract only the transverse part of the full vertex. In order to project the quantity above 
on the tree-level vertex, the following function will be evaluated 

(P ' *' ^ «^(p, q, k) D«i(p) D^q) D%{k) l£^(p, q,k) M 

= T% L p A3 > abc (p,q,k) Gf^ b %p,q,k) 

j«Lfi,A»,abcfa g) k) p M p va{q) p pUk) ptl ; L,A» 1 .6c (pj g> fc) D(p] D{q) D[k) > 



G A „ p abc {p, q,k) = -< Al(p)AUq)A p (k) > . (18) 



where the D's are the scalar gluon propagators. Clearly, 
this function is equal to 1 if the full vertex is equal to 
the bare vertex. Note that, by contracting the Green's 
function with the lattice version of the tree-level vertex 
instead of the continuum version corrects for (relevant) 
discretization effects. Also note that the exponential pre- 
factor in Eq. JSJ implies a pre-factor exp[— iir(P( Jl +p\)/N] 
for the gluon propagator D^f(p). Thus, considering Eq. 
(|15fl . all these pre- factors cancel each other both in the 
numerator and in the denominator of Eq. I|20[) above. 

The normalization factor in the denominator is in 
some kinematical cases quite simple, but in general very 
lengthy and will not be given here explicitly. Also, this 



normalization factor vanishes for the largest momentum 
in each momentum configuration due to the cosine fac- 
tors. Thus, this momentum cannot be considered. The 
same applies to the case where all momenta vanish. 

Let us stress that, although the quantity in Eq. i|21|) 
is 1 if the full and the tree-level vertex coincide, it will 
in general have contributions from tensor structures not 
appearing at the tree level. The situation would be even 
more complicated in the case of a vertex not totally anti- 
symmetric in color space. In the continuum, the general 
tensor structure of the total color-antisymmetric part of 
the full three-gluon vertex is given by |53| 



r fvp abC (P, 2, k ) = -if abC {MP, <1, k)S^(p p - q p ) + B(p, q, k)S^{p p + q p ) + C(p 7 q, k){p u q p - S^p ■ q)(p p - q p ) 

+ S (ppq^ku +Puq P k p )) +cyc. perm. , (22) 
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where A, B, C and S are scalar functions. (See again Ref. y^j for the symmetry properties of these functions.) Then, 
one can verify that, in the continuum, the quantity l|21[) contains only contributions from the functions A and C. For 
example, in the orthogonal configuration, using the above equation one finds 

G A " C (p, q, |) = {xyA(p, q, k) + (x + y) [2yA(k,p, q) + 2xA{q, k,p) 

+ (xy + 2y 2 )C(k, Pl q) + (2x 2 + xy)C(k, q,p) + xyC(q,p, k)]}/{2x 2 + bxy + 2y 2 ) , (23) 



where x = q 2 and y = p 2 . The expression is lengthier for 
other kinematical configurations. 

It then only remains to determine the 3-point Green's 
function. This can be done in a straightforward way 
using the definition (|18fl an d the Fourier-transformed 
gauge fields defined in Eq. ©. The only problem is 
that the functions A a ^(p) vanish on average. This induces 
very large fluctuations in the calculation of the Green's 
function, particularly at large momenta. This problem 
makes the extraction of the equal-momentum configu- 
ration quite complicated. In the orthogonal configura- 
tion the situation is slightly better, since some field com- 
ponents identically vanish even for non-vanishing mo- 
menta, as discussed in Section IIIII reducing the statis- 
tical noise. Also note that, in the evaluation of Eq. (|21|l . 
the Green's function and the propagators are calculated 
independently and divided after averaging over all lat- 
tice configurations. The error can then be determined by 
error propagation, but due to the smallness of the statis- 
tical errors on the propagators, their errors are neglected 
here. 

Let us note that, from the discussion in Section 
IIII Al we found that (with our notation) the contin- 
uum momentum-space gluon field has mass dimension 
— 1 — d/2 while the gluon propagator, again in momen- 
tum space, has mass dimension —2. This implies that 
the full Green's function defined in Eq. i|18fl has mass di- 
mension —3 — d/2. The same result can be obtained by 
considering Eq. (|19|l after observing that the three-gluon 

vertex ^\ a '^ e ^ {p, q, k) has mass dimension 3 — d/2 [see 
also the tree-level result in Eq. (|14|) ]. Thus, the quantity 
considered in Eq. I|2(J|) is inherently dimensionless and 
we do not have to multiply it by any power of the lattice 
spacing a. On the other hand, in order to get the corre- 
sponding continuum quantities we have to multiply the 
momentum-space lattice gluon field and the gluon prop- 
agator by \fj3 and by (3, respectively. It follows that the 
scalar function G A (p, q, <fi) has to be divided by /3 3 / 2 in 
order to obtain the corresponding continuum quantity. 

Three different momentum configurations are shown 
in Fig. E| The first is the momentum configuration used 
in 123 . 130L l3l| to study the behavior of the strong- 
coupling constant. Note that this configuration is not 
well-defined in the continuum limit in case the Gribov- 
Zwanziger scenario is correct, i.e. if the gluon propagator 
vanishes at zero momentum. Of course, on a finite lattice 
this is not a problem. In the second case, two momenta 
are equal and orthogonal. Finally, the third case cor- 



responds to three momenta with the same magnitude. 
Furthermore, the behavior of G A (p, q, 4>) as a function 
of p and q is shown in Fig. EI 

As said above, the statistical errors are quite large 
at large momenta, even after averaging over more than 
20, 000 configurations. Let us recall that, in some cases, 
the data shown in Fig. [S] have been obtained by sum- 
ming more than 700 different terms. In particular we 
checked that, at large p, the central values of the data can 
change substantially if one considers different subsets of 
the whole set of configurations used in the analysis. This 
suggests that the statistical error obtained at large p is 
likely under- estimated. 

Despite the large statistical fluctuations at large mo- 
menta, the vertex clearly decreases in the limit of small 
momenta. This behavior can be observed in the three 
plots in Fig. [H] and also in Fig. for all directions ap- 
proaching the limit p = q = k = 0. On the other hand, 
at the smallest momentum point on the 30 3 lattice at 
(3 = 4.2 (see top panel in Fig.|3Jl, i.e. for the orthogonal 
configuration, the vertex is clearly negative. Of course, 
with our data we cannot say if the vertex would stay fi- 
nite or would become larger (in absolute value) as the 
zero momentum limit is approached. Indeed, we know 
that, with our lattice volumes and (3 values, the true 
infrared regime is not reached yet, since the gluon prop- 
agator is not suppressed at small momenta ■ Clearly, 
if the 3-point Green's function stays constant in the in- 
frared limit, while the propagators get suppressed, the 
vertex F A v '// bc (p,q,(j)) may be enhanced, as predicte d by 
studies using functional methods in four dimensions [TfJ 
and in three dimensions |ll| . 

Finally, in order to test whether other tensor structures 
are relevant, the vertex has also been contracted with 
itself instead of using the tree-level vertex. Results are 
not presented here, but we find that the corresponding 
scalar function shows a strong increase compared to the 
tree-level behavior at large momenta. This suggests that 
other tensor structures also contribute to the vertex in a 
non-trivial way. 



B. Ghost-gluon vertex 

The ghost-gluon vertex can be treated essentially along 
the same lines as the three-gluon vertex. Let us recall 
that at the tree level, in the continuum, this vertex is 



Three-gluon vertex, all momenta equal 




FIG. 8: The scalar function G A c (p,q,(/>) defined in Eq. 
(12( )ti _ Full symbols correspond to (5 = 4.2 and open sym- 
bols to P = 6.0; circles are used for V = 20 3 and triangles for 
V = 30 3 . In the top panel we show results for the orthogonal 
configuration with one momentum (q) vanishing. In the mid- 
dle panel we consider an orthogonal configuration with two 
momenta having the same magnitude (p = q). In the bottom 
panel we plot data for the case with the three momenta equal 
(p = q = k). 



FIG. 9: The scalar function G A c (p, q, 0) defined in Eq. 
as a function of the magnitude of the gluon momenta p and q 
(for the orthogonal configuration). Here, for each data point, 
we plot only the central value. The data for (3 = 4.2 and 
P — 6.0 are plotted together and they are interpolate d by 
Gouraud shading, as implemented in the ROOT package [54j|. 
In the top figure we used data for the lattice volume V = 20 3 ; 
in the bottom one we consider the lattice volume V = 30 3 . 
Spikes indicate positions where due to fluctuations the value 
is outside the drawing range. The spike at p — q = is 
an artifact since this quantity cannot be evaluated when the 
three momenta are null. 



given by 

rf™ A < abc (p,q,k)=igf abc qfi , (24) 
while on the lattice one finds 

r tl,L,ccA,abcfo fc) = lgf abc el ,kJN^ CQs( ^ 

Again, the cosine and the exponential are lattice arti- 
facts, going to 1 in the formal continuum limit a —> 0. 
Also, on the lattice it is only possible to determine the 
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full Green's function 

Gf A ' abc ( P , g, fc) = i < c a (p)c b (q)A^k) >, (26) 
I 



where D G is the ghost propagator and we used Eq. (|25Jl . 
Note that we have neglected in the denominator the color 
off-diagonal components of the propagators and used the 
relation f abc f abc = N C {N 2 ~ 1). Also, due to the im- 
plicit contraction of the full Green's function with a gluon 
propagator, only one tensor structure survives, which is 
proportional to the incoming anti-ghost momentum q. 
Finally, as for the three-gluon vertex, the exponential 
pre-factors cancel out in the numerator and in the de- 
nominator of the scalar function defined above. 

For the ghost-gluon vertex, the normalization factors 
are quite simple and will be given explicitly. In the 
orthogonal case, we considered a gluon momentum k 
aligned along the y-axis, while the incoming anti-ghost 
momenta q is chosen along the x-axis. This implies 
qfiPfiv(k) = q v . Thus, the denominator in Eq. I|29[) is 
proportional to 

$ cos2 (^) = fx cos 2 ^), (30) 

A" 

where the last equality follows because the momentum q 
is aligned along the x-direction. Clearly, in this case one 
recovers the kinematical normalization considered in Rcf. 
|25| . In the case of equal momenta, the situation is more 
complicated and an explicit Gram determinant appears. 
Indeed, by choosing the gluon momentum k in the x — y 
plane and the incoming anti-ghost momentum q in the 
x — z plane, the denominator is proportional to 

ql COS 2 (g x ) (\ - \ + C0S \q z )ql . (31) 

In the case of equal components q x — q z — k x = k y = q 
this expression simplifies to 3q 2 cos 2 (q)/2. Note that, in 
both cases, the normalization would be the same in the 
four-dimensional case, since all momenta considered are 
in a three-dimensional sub-space. Finally, it should be 
noted that the only Lorentz indices of G^ A ' abc (p,q,k) 
contributing to the numerator are those corresponding 



where c a (respectively c b ) is the ghost (respectively anti- 
ghost) field. Then, the scalar quantity we evaluate is 
defined by 



(27) 



(28) 



(29) 



to non-zero components of the incoming anti-ghost mo- 
mentum q. Thus, in the orthogonal configuration the 
argument of the exponential pre-factor appearing in Eq. 
(|5|) is always zero for the Lorentz components contribut- 
ing to Gf A ' abc {p,q,k). 

Let us also note that, when contracted with the trans- 
verse projector P M! ,(fc) of the gluon momentum k (as done 
here) , the vertex should be invariant under the exchange 
of the ghost and of the anti-ghost fields |2J|. This is of 
course true at the tree- level [see Eq. (EU], since it corre- 
sponds to an exchange of the color indices a and b and to 
replacing q by p — — k — q. Indeed, the term proportional 
to fcp vanishes, due to the contraction with P^ u (k), and 
the minus sign of — q^ cancels with the minus sign related 
to the antisymmetry of the structure constant. In the 
general case this invariance is a consequence of a global 
SL(2,R) symmetry between ghosts and anti-ghosts (see 
Appendix A in Ref. [|J for a proof of this symmetry in 
Landau gauge). As the Green's function (|33|l is implic- 
itly contracted with a gluon propagator, it follows that 



Gf A ' abc (p, q, k) = Gf A - bac (q,p, k). (32) 

This result is also related to the symmetry of the 
Faddeev-Popov matrix under simultaneous exchange of 
space and color indices. 

The only remaining task is then to determine the quan- 
tity defined in Eq. (|26[) above. Since on the lattice 
one does not consider ghost (and anti-ghost) fields, the 
(equivalent) expression to be used is |25| 

Gf A > abc (p, q,k) = ±< (M-Y'iP, 9) ^{k) > . (33) 

One can employ the point-source method, used for evalu- 
ating the ghost propagator, also in the case of the ghost- 
gluon-vertex function. Indeed, by using translational in- 
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rf L < ceA > abc (p, q, k) Gf A > abc (p, q, k) 



itl, L,ccA t abc 



(P,q,k) D- G d {p) D b £{q) D%(k)T 



-\t\,L,ccA.def 



(p, q, &) 



T ti ' L ' c5A ' abc (p,q,k) G 



ccA,abc 



(p, q, k) 



Tf L ' cSA ' abc {p,q,k)P flv (k) TV 



t\,L,ccA,abc 



(p,q,k) D G {p) D G {q) D(k) 



-i<7/ ab % cosfe) Gf A > abc (p,q,k) 



g 2 N c (N? - l^cos^P^fc^cos^) D G {p) D G {q) D(k) ' 

I 
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Ghost-gluon vertex, one momentum vanishing (PS) \ Ghost-gluon vertex, one momentum vanishing (PWS) \ 




5?" 



QA A AqA, 



3.5 
q [GeV] 



I I I I l_l I I I I I l_l I I I I L_l I I I I I l_l I I I I l_l I I L 



0.5 



1.5 



2.5 



3.5 
q [GeV] 



Ghost-gluon vertex, orthogonal momenta with two equal (PS) I Ghost-gluon vertex, orthogonal momenta with two equal (PWS) 



o- 



2.5- 
2 

.1.5 



0.5- 
- 



f 



I . i . . I 



I ... i I 



57 



2.5- 

2 : 
.1.5 



0.5- 
- 



if ,L .1°* ^ ^^'^OA O AgA i 



0.5 1 1.5 2 2.5 3 3.5 

q [GeV] 



J i i_i i I i i i i_l i i i i l_i i i i I i i_i i L 



0.5 1 1.5 2 2.5 



J I l_l I I L 



3 3.5 

q [GeV] 



Ghost-gluon vertex, all momenta equal (PS) 



Ghost-gluon vertex, all momenta equal (PWS) 




'a a ao 



'a 6£ 



■o.s. r i i i 



_i i i i i i i i i i i i i i_ 



5 

q [GeV] 



FIG. 10: The scalar quantity G cSA (g,fc,0) defined in Eq. I^TJ as a function of the magnitude of the incoming anti-ghost 
momentum q. Full symbols correspond to (3 — 4.2 and open symbols to f3 — 6.0; circles are used for V = 20 3 and triangles for 
V — 30 3 . Results on the left-hand side have been obtained using a point source (PS) while on the right-hand side we considered 
a plane- wave source (PWS). In the top panel we show results for the orthogonal configuration with the gluon momentum k 
vanishing. In the middle panel we consider an orthogonal configuration with the two momenta (q and k) having the same 
magnitude. In the bottom panel we plot data for the case with the three momenta equal. For the number of configurations 
considered in the PWS case, see caption of Table ITTT1 
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I Ghost-qluon vertex, orthogonal momenta, V=20 3 




Ghost-gluon vertex, orthogonal momenta, V=30 3 | 




FIG. 11: The scalar quantity G c5A {q,k,(j>) defined in Eq. 
1271 as a function of the magnitude of the incoming anti-ghost 
momentum q and of the gluon momentum k (for the orthog- 
onal configuration) using the point-source method. Here, for 
each data point, we plot only the central value. The data for 
j3 — 4.2 and (5 = 6.0 are plotted together and they are inter- 
polated by Gouraud shading, as implemented in the ROOT 
package |54|. In the top figure we used data for the lattice 
volume V = 20 3 ; in the bottom one we consider the lattice 
volume V — 30 3 . Spikes indicate positions where, due to fluc- 
tuations, the value is outside the drawing range. The large 
spikes below the smallest non-zero anti-ghost momentum q 
are an artifact of the interpolation. 



variance we can write 8 

< M' 1 (x, y) A(z) > = < M _1 (0, y - x) A{z - x) > . 

(34) 



Here, in order to simplify the notation, we omit the color indices. 



Then, by using momentum conservation, i.e. p = —k — q, 
and the equation above we obtain 

I^ e 2-(-(fe+?)x+ TO +M/A' < M~ x (x,y)A(z) > 

xyz 

— 1_ e 2irt[q{y-x)+k(z-x))]/N 
xyz 

x < M _1 (0, y - x)A{z -x)> 
= Y^ e 2ni{qv+kz)/N < M- 1 {Q,y)A{z) > 
= J2e 2m(qy+kz)/N < M- 1 (x,y)A(z)S x0 > 

xyz 

= ^2 e 2ni( qy+k *)/N < M - 1 (x,y)A(z)(S x0 - I) > . 

xyz 

In the third line we have re-defined the indices y and z 
and summed over x. In the last line we added a term that 
vanishes for non- vanishing ghost momenta |46j . Clearly, 
the quantity S x q — 1/V has a null value when summed 
over x, i.e. the inversion of the Faddeev-Popov opera- 
tor is done in the subspacc orthogonal to the (trivial) 
kernel of M. Thus, one can evaluate the Green's func- 
tion by considering the Fourier-transformed gluon field 
and by inverting the Faddeev-Popov operator using the 
point-source method. The final step requires one to eval- 
uate the Fourier transform using the incoming anti-ghost 
momenta q, i.e. the one appearing in the tree-level ver- 
tex. The third momentum (i.e. p) is implicitly defined 
by momentum conservation. As said above, the point- 
source method has the advantage of only N% — 1 inver- 
sions per configuration but, compared to the plane- wave- 
source method, requires many more configurations in or- 
der to achieve a given statistical accuracy. 

Note that, since the momentum-space ghost propaga- 
tor Dc has mass dimension —2, we have that the full 
Green's function G^ A ' abc (p,q,k) [see Eq. (|33|l ] has mass 
dimension — 3 — d/2. At the same time, the ghost-gluon 
vertex r* 1 ' L * CCj4 ' o6c (p, q, k) [see Eq. (j23J] has mass dimen- 
sion 3 — d/2. Thus, also in this case the scalar function 
considered is clearly dimensionless and we do not have to 
multiply it by any power of the lattice spacing a. On the 
other hand, in order to get the corresponding continuum 
quantities we have to divide it by /3 1 / 2 . 

Results are shown in Figs.^Jand^J I n the left column 
of Fig. we show the results obtained using the point- 
source method, while on the right column we present the 
data obtained using the plane-wave source, for the ver- 
tex and for the ghost propagator appearing in Eq. 129|) . 
We see agreement between the two sets of data but, as 
said above, the fluctuations are much smaller in the lat- 
ter case, the relative error being of the order of 5% for 
all momenta. On the other hand, when considering the 
point-source method (left column) the fluctuations be- 
come rather large with increasing gluon momentum, i.e. 
the results rapidly lose accuracy. Let us finally note that 
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I Ghost-qluon vertex, orthogonal momenta, V=20 3 




Ghost-gluon vertex, orthogonal momenta, V=30 3 | 




FIG. 12: Same as Figure 1111 but using a mixed point- 
source/plane- wave-source method. The number of configu- 
rations considered here are 283 and 269 for the lattice volume 
V = 20 3 at (3 = 4.2 and (3 = 6.0, respectively, and 227 and 
213 configurations for V = 30 3 at (3 = 4.2 and (3 = 6.0. 



the results on the right column do not change visibly if 
one evaluates the ghost propagator using the point-source 
method and the vertex using the plane-wave source (see 
Figure EJ . 

Clearly, in all cases the vertex is essentially constant 
and of order one, for all momentum configurations, con- 
firming the results obtained in the 4d case |25J, 26] . This 
would imply a fulfillment of the corresponding Slavnov- 
Taylor identity (5{|. This result is also in very good 
agreement with the functional predictions in the 3d case 
|llL Il2j | and the usual assumptions made in functional 
calculations. 



V. SUMMARY AND OUTLOOK 

In this work we have evaluated the three-point vertices 
of (pure) Yang- Mills theory using a wide range of momen- 
tum configurations. In order to reduce the computational 
cost we considered the 3d SU (2) case. We note that the 
data for the three-gluon vertex suffer from large statisti- 
cal errors at large momenta, even when considering more 
than 20,000 configurations. Let us mention that, from a 
preliminary study |56j, this problem seems to be less se- 
vere in the 4d case. For the ghost-gluon vertex the data 
also show a large ratio noise/signal if one uses the point- 
source method. In this case the use of the plane-wave 
source (for the vertex) allows to reduce the error below 
10% already with a few tens of configurations. 

As for the infrared behavior of the vertices, we found 
that the three-gluon vertex becomes very small as the 
momentum decreases. This mid-momentum suppression 
is similar to the behavior considered for this vertex in 
Ref. [j| in order to obtain a positive semi-definite gluon 
propagator. From our data it is difficult to say what 
would be the behavior of this vertex at very small mo- 
menta. However, in one of the kinematical configurations 
considered here, the three-gluon vertex becomes negative 
at the smallest non-zero momentum (about 240 MeV). 
Thus, a possible scenario could be a vertex becoming 
larger (in absolute value) as the zero momentum limit is 
approached. Let us note that a positive infrared diver- 
gent three-gluon vertex has been recently found in func- 
tional studies in three dimensions |ll| . 

The ghost-gluon vertex stays constant and essentially 
equal to the tree-level value in the range of momenta 
considered. This is agreement with various theoretical 
predictions llOL IllL [l2j and with numerical results for 
the 4d case [H, |26| . 

Thus, our results, albeit exploratory, seem to support 
(at least at the qualitative level) the Gribov-Zwanzigcr 
and Kugo-Ojima scenarios of confinement and the cen- 
tral assumptions usually considered in functional meth- 
ods. In particular, we have shown (see Figure UJ that 
the smallest non-zero eigenvalue of the Faddeev-Popov 
matrix goes to zero in the continuum limit, i.e. a (con- 
tinuum) Landau configuration should belong to the first 
Gribov horizon. On the other hand, these results should 
be taken with caution, as we know that in 3d with the 
(physical) volumes considered here the true asymptotic 
infrared region has not been reached yet [l^ . 

We plan to extend this study to larger lattice volumes, 
in order to explore the far infrared limit. We also plan 
to consider other gauge groups, especially the physical 
SU(3) group. Finally, one should consider possible sys- 
tematic effects related to the existence of Gribov copies. 
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